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ABSTRACT 


A simple mathematical model was developed to simulate 
the closed portion of the cycle for a quiescent chamber 
compression ignition engine utilizing the assumption of 
perfect gases and the first law of thermodynamics. Various 
input paramentes were used in trend analysis to check the 
model. The output from the computer program was Compared to 
test data from a four inch bore, open chamber semi-quiescent 
diesel engine run at the Sloan Test Laboratory. This 
computer model was then modified to simulate the expansion 
stroke of a newly developed, two cylinder cycle 
reciprocating engine. The model was then run to determine 
the optimum point of fuel injection for the new engine. 
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Chapter |] 
INTRODUCTION 


The compression ignition (CI) engine is thriving in new 
found popularity amongst automobiles, medium-duty and 
heavy-duty freight transport trucks, marine propulsion and 
auxiliary systems, and various other industrial 
applications. The United States Department of Energy (DOE) 
has recently predicted that diesel fuel consumption will 
exceed gasoline consumption in this country by the year 
2000. This is primarily due to the shift to diesel power in 
the automotive and truck freight transport industries to 
take advantage of the high efficiency, high power-to-weight 
ratio engines. The present daily diesel fuel consumption of 
the order of 108 liters is expected to increase by as 
much as 50 percent by the turn of the century. The 
increasing fuel consumption rate, coupled with the constant 
concern of diminishing oil reserves, has prompted renewed 
interest in improving the operating efficiency of the 
conventional compression ignition engine. Any small 
improvement in engine efficiency will obviously result in an 


enormous savings in petroleum. 





The approaches currently pursued to improve compression 
ignition engine efficiency include increasing the 
compression ratio and the development of the "adiabatic" 
engine. The former involves turbocharging and the 
improvement to piston ring technology. The latter approach 
concentrates the most emphasis on insulating the engine. 
This requires the use of temperature-resistant ceramic 
cylinder liners for combustion cylinders whose gas wall 
temperatures can be of the order ot 1200 degrees Kelvin. In 
addition to these two approaches, there are many other 
avenues of research in progress that involve improvements 
that will increase the compression ratio, decrease the heat 
loss from the engine, or increase the combustion efficiency 


through improved combustion chamber design. 


Instead of improving upon the conventional compression 
ignition engine, a new cycle engine design is under 
development. This new design, proposed by Carmichael 2’, 
consists of a two cylinder cycle which divides the functions 
of a conventional four-stroke diesel cycle into two parts. 
The new engine has one cylinder which compresses the 
incoming air charge and another cylinder which acts as the 
combustion chamber and expansion cylinder. These two 
cylinders are interconnected by a regenerative heat 
exchanger. The regenerator acts as the heart of the new 
design. Through the use of new ceramic materials, the 


regenerator will act as a heat transfer medium by 
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transfering a portion of the heat from the exhaust gases to 


the incoming air charge. The temperature of the incoming 
air will be elevated twice. The first temperature increase 
is due to the compression process in the first cylinder. 
This cylinder, in turn, will transfer its air charge through 
the ceramic matrix of the regenerative heat exchanger, thus 
boosting the temperature for the second time. After passing 
through the regenerator, the incoming air charge will be of 
sufficient temperature to accomodate spontaneous combustion. 
with this high temperature, the high compression ratio of 
the conventional compression ignition engine is not required 
to obtain work from the cycle. Figures ] and 2 depict the 
pressure vs. volume and temperature vs. volume diagrams for 
the new cycle as compared to a conventional diesel cycle. 
The preliminary design of this new two cylinder cycle engine 
indicates that an improvement to thermal efficiency can be 


achieved over the conventional compression ignition engine. 


An important element of the engine design process is 
the capability to predict, with an acceptable degree of 
accuracy, the energy release during combustion as a function 
of time. This process is extremely complex in that it 
involves the injection and atomization of fuel, the 
evaporation and mixing of the fuel with the air charge, 
foliowed by the various phases of combustion. The ability 
to accurately predict the heat release rate is vital to the 


engine designer when analyzing a new engine design. 
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This thesis is an attempt to assimilate various diesel 
engine combustion models to produce a simple, yet accurate, 
model to be used in the continuing evaluation of the new two 
cylinder cycle reciprocating engine. The proposed model can 
be utilized on a personal computer system to determine Che 
optimum point of fuel injection for the new engine. The 
model has the capability to evaluate two different fuel 


types (i e. iso-octane and propane). 
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Chapter 2 
DESCRIPTION OF THE TWO CYLINDER CYCLE 


RECIPROCATING ENGINE 


The Two Cylinder Cycle Reciprocating Engine consists of 
one compression cylinder and a pair of combustion/expansion 
cylinders, see figures 3 and 4. The compression cylinder 
consists of an intake valve and two exhaust valves, One to 
each expansion cylinder. Each expansion cylinder has its 
own fuel injector. The regenerator cavity contains an 
exhaust valve in addition to the ceramic matrix regenerator. 
The pistons of both the compression cylinder and the 
expansion cylinders are considered to be of simple geometry 
with flat heads.» The expansion cylinder piston incorporates 
no unique features to increase turbulence or swirl, thus it 
is similar to a direct injection, quiescent chamber diesel 
engine cylinder. The five valves and three pistons are 
actuated by a camshaft that allows the compression piston to 
Operate at twice the speed of an expansion piston. The 
compression cylinder will alternately provide a compressed 
air charge to each expansion cylinder via the regenerator. 
A typical cycle can be illustrated by referring to figures 3 


and 4. 
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has just closed (it only opened for a very short time just 


before the compression cylinder reached TDC), valve 3A just 
opened, valves 2A and 3B are closed. (This is the same as 


step l except that cylinders À and B are reversed.) 


Step 6: The compression cylinder is approximately 90° 
ATDC and is in the middle of an air charge induction. 
Cylinder A is still exhausting and cylinder B is in the 
expansion process. Valve l is open, valves 2B, 3B, 2A are 
closed and valve 3A is open. (This is the same as step 2 


except that cylinders A and B are reversed.) 


Step 7: The compression cylinder has just reached BDC 
and has completed induction of an air charge. Cylinder A is 
exhausting through the regenerator and cylinder B is still 
expanding. Valve 1 just closed, valves 2B, 3B, 2A are 
closed and valve 3A iS open. (This is the same as Step 3 


except that cylinders A and B are reversed.) 


Step 8: The compression stroke is approximately 90% 
BTDC and in the middle of compressing an air charge. 
Cylinder A is completing its exhaust stroke and cylinder B 
is ending its expansion stroke. Valve 1 is closed, valve 
2B, 3B and 2A are closed, valve 3A is still open. (This is 
the same as step 4 except that cylinders A and B are 


reversed.) 


aA. 





Step 9: This is the same as step 1. 


Figures 5 and 6 show the temperature and pressure as a 


function of cylinder volume for a cycle. 


The table below summarizes the sequencing of the valves 


for a complete cycle of an expansion cylinder. 


Tab : uenci alves 
Valve 

1 2A 3A 2B 3B 
Step 1: X X X X 0 
Step 2: O X X X O 
Step 3: X X X X 0 
Step A: X X X X O 
Step 5: X X O X X 
Step 6: O X O X X 
Step 7: X X O X X 
Step 8: X X O X X 
Step 9: X X X X O 


where X = Valve closed 


and O Valve open 


As can be readily seen, the valve timing sequence is 
rather complex. The timing sequence must be such as Co 


allow the impulse transfer of the air charge to occur 


15 





without possible blow-down to the atmosphere or charging the 
wrong cylinder. A shift of the crank angle must be 
considered to optimize the air charge transfer sequence to 
the on-line expansion cylinder. Thus, the valve timing 
sequence is a critical factor in the correct and efficient 
Operation of this new engine design and must be dealt with 


appropriately. 
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Chapter 3 
COMBUSTION AND COMBUSTION MODELING 


(An overview) 


3.1 Description of Diesel Engine Combustion 


The diesel engine combustion process is exceedingly 
complex and not very well understood. Combustion in the 
diesel engine is characterized by compression ignition, a 
non-uniform fuel and air distribution in the combustion 
chamber, and a continuous mixing throughout the period in 
which combustion occurs. Due to the initial conditions in 
the chamber when fuel is first injected, the air charge in 
the cylinder is of sufficient temperature and pressure to 
support a chain-reaction. However, combustion in the 
compression ignition engine is governed by the local 
conditions in each part of the charge and not dependent on 
the ودک‎ of the flame from one point to another. 
Therefore, the rate of combustion is dependent on the state 
and distribution of the fuel and upon the pressure and 


temperature within the cylinder. ?! 


3.1.1 The Phases of Combustion 


Sa 





Ricardo described the diesel engine combustion process 
as taking place in three stages; namely the delay period, a 
period of rapid combustion, followed by burning at a 


(3) Lyn ii described the burning process in 


controlled rate. 
three slightly different phases. The first phase is a 
period of rapid combustion which lasts for only three 
degrees crank angle. The second stage is characterized by a 
decreased rate of heat release lasting approximately 40 
degrees crank angle. The third period consists of the fuel 


burning at a very slow rate which may persist through the 


remainder of the expansion stroke. 


A combination of the descriptions of Ricardo and Lyn 
may be more appropriate. The stages of combustion could be 
divided into ignition delay, premixed burning, diffusion 
controlled combustion and the tail of combustion. (?'6! 


Figure 7 depicts the four stages of combustion in a heat 


release diagram. 


3.1.1.1 Ignition Delay 


The term ignition delay, or ignition lag, describes the 
time required by the preliminary reactions that occur prior 
to the appearance of flame. The ignition delay is broken 
down into a physical delay and a chemical delay. The 
physical delay period occurs between the beginning of fuel 
injection and the onset of chemical reactions. During this 


period, the fuel is atomized, vaporized, mixed with air and 


zd 





raised in temperature. This process is sometimes 
collectively referred to as preparation. The chemical 
delay period immediately follows the physical delay period 
and terminates at inflammation or ignition. This period is 
characterized by chemical reactions starting slowly with 


pre-flame oxidation of the fuel followed by local ignition. 


The ignition delay will vary according to cylinder 
temperature, cylinder pressure, the type of fuel, the 
initial temperature of the fuel, the characteristics of the 
fuel injectors and the turbulence in the cylinder. The 
physical delay is small for light fuels but can become the 
controlling factor for heavy, viscous fuels. The physical 
delay can be significantly reduced by using high injection 
pressures and high turbulence Eo expedite the breakup of the 


fuel jet. 


Semi-empirical relationships have been developed to 
describe the ignition delay. An estimate for igniton delay 
was developed by Wolfer in 1938: 17? 

t = 0.44P7 +19 exp (4650/7) 
where: t = ignition delay in milliseconds 
P = cylinder pressure in atmospheres 


and T = temperature in degrees K at 


ignition. 
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An estimate by Clarke! in 1970 is quite similar to that by 
Wolfer: 
t » 0.22exp(5500/7) PO - 777 
where: t = ignition delay in seconds 


i 


cylinder temperature in degrees K 


and P 


cylinder pressure in Nm. 


Still another empirical expression for ignition delay was developed 
by Spadaccini and Tevelde*?? from experiments for NASA in 1979 with 
diesel fuel in a steady flow facility: 

t = 2.43x10P”Zexp(41560/RT) 
where t = ignition delay in seconds 


P 


pressure in atmospheres 
T = mixture temperature in degrees K 
and R = gas constant in atm cm?/gmole°K. 
Figure 8 represents the effects of temperature and pressure on ignition 


delay as determined from the estimates by Wolfer. The Spadaccini and 


Tevelde and Clarke relationships yield somewhat similar results. 


When uSing ignition delay expressions, it must be emphasized that 
differences in engines, fuel properties (especially cetane number), fuel 
injectors and actual engine temperatures and pressures make the 
calculation rather approximate. These formulas are also very limited by 
their use of bulk temperatures, with no consideration of local 


compositions or temperatures, (10? 


EU 





3.1.1.2 Premixed Burning 


In the premixed burning stage, flame occurs at one or more 
locations and spreads turbulently. The rate and amount of combustion 
during this stage is directly related to the fuel preparation rate and 
the length of the ignition delay period. Since this stage of combustion 
is one of premixed combustion, little carbon (soot) is produced 
resulting in little radiation heat transfer. However, since the 
combustion rate is so intense, combustion generated noise is controlled 


(11) 


by this stage of combustion. Figure 9 depicts premixed burning in a 


cylinder. 


3.1.1.3 Diffusion Controlled Burning 


Once the prepared, or premixed, fuel has burned, the combustion 
process slows down. The combustion rate in this stage will be dominated 
by the rate of local air entrainment. Since the temperature in the 
cylinder is favorable for ignition in this stage, the air/fuel mixing 
process will control the rate of combustion. This preparation of the 
fuel will be governed by the turbulence and swirl in the cylinder. 
Lyn 4? estimated that approximately 40 percent of the heat release from 


the combustion of fuel comes from this stage. Figure 10 shows the 


diffusion burning process in a cylinder. 


3.1.1.4 Combustion Tail Stage 


This last stage of combustion is characterized by the cylinder 
pressure and temperature falling as the expansion process continues. 
The rate of combustion tails off due to the chemical kinetic effects as 
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the chemical reaction rate slows. In this stage, the reaction rate will 
become the controlling factor instead of the air/fuel mixing process. 
This stage is also characterized by diffusion combustion with a high 
production and combustion of soot particles with a resultant high rate 
of radiation heat transfer. This last stage of combustion can proceed 
through the completion of the expansion stroke and can contribute upto 


{4} 


20 percent of the total heat release. Figure 11 represents a typical 


heat release rate diagram showing the four stages of combustion. 


3.2 Combustion Modeling 


The combustion process is often considered the most important 
aspect of an internal combustion engine, but, at the same time, the 
least understood and most complex. A mathematical model depicting 
combustion would require good models of the fuel system to include the 
injection/fuel pump, the injector nozzles, and fuel lines. 
Additionally, models of fuel atomization, vaporization, fuel/air mixing, 
cylinder air motion, chemical kinetics and pre-mixed and diffusion 
mixing would be required. A model as comprehensive as this has yet to 


(12) states that this type of "combustion 


be deveoped. Spaulding 
modeling is impossible." He justifies this by pointing out that the 
number of governing restraints and rules outnumber the degrees of 
freedom and, in addition, the requirements of low cost, speed and 
accuracy must also be met. Since the complexity of the real combustion 


process is so overwhelming, substantial simplifying assumptions must be 


made to obtain solutions. 
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3.2.1 Types of Models and Uses 


{13} categorized combustion models into three categories 


Bracco 
based on their uses in examining different engine problems. The 
categories are the zero-dimensional (or thermodynamic) model, the 
quasi-dimensional (or entrainment) model, and the multi-dimensional (or 


detailed) model- 


3.2.1.1 Zero-dimensional Mode 


The zero-dimensional model is structured around a thermodynamic 
analysis of the engine cylinder contents during the cycle. The 
assumptions include one-dimensional flow, isentropic adiabatic flow 
through nozzles simulating flow past valves, and unburned mixtures as 
mixtures of air, fuel vapor and residual gases. Specific heats of the 
gas mixture are modeled using polynomial functions of temperature. 
Compression is assumed to be adiabatic. Combustion assumes 
thermochemical equilibrium and progressive burning via mass elements, 


The expansion process assumes thermochemical equilibrium. 


Heat transfer is modeled using correlations between the Nusselt, 
Prandtl- and Reynolds numbers from heat transfer in steady turbulent 
flow over flat plates and pipes. These relationships are in the form 
Qu 

Nu = aRePpr° 
where a, b, and c are obtained from experimental data for a specific 


engine. 
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The combustion process is generally modeled from an apparent heat 
release or an experimentally obtained fuel burning rate. One of the 
most widely used correlations is based on the Wiebe Function. In this 
function, the fuel burned is expressed as a fraction of the total fuel 
injected, ‘>! 

FB = 1 - expl-K (t) F1 
where FB = fracticn of fuel burned/total 
injected 
t = time from ignition 
K; = shape factor for combustion curve 


K, = combustion efficiency coefficient. 


Another typical function form is the cosine function: tH? 
X(@) = (1/2)11 = cos r [(© = GC E 
where X(8) = mass fraction burned at crank 
angle ٥ 
e = crank angle at the start of 


combustion 


and ^6, z burn duration. 


There are numerous other combustion models that utilize various 
heat release patterns. Some replace the heat release curve with two 
Straight lines. In this type of combustion model, one line simulates 
the rapid combustion of the bulk of the injected fuel and the other line 
represents the slower combustion phase further down the expansion 


Stroke. 


OA 
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An empirical model developed by Whitehouse and Way 14! is based on 
elementary combustion principles. Fuel is assumed to be prepared for 
combustion as a result of fuel-air mixing. The reaction rate calculates 
the burn rate in the premixed stage of combustion. The preparation rate 
becomes governing during the diffusion burning phase as the fuel is 
assumed to burn as rapidly as it is prepared. (The Whitehouse and Way 


model will be dealt with in detail in a later chapter.) 


In general, thermodynamic combustion models are useful when 
performing a design trade off or comparison analysis to evaluate the 
effects of change in engine design and operation. Since, however, the 
details of the combustion process are an input to the model, the results 
can only indicate what will transpire if the engine burns in the 
specified manner. These models cannot address the feasibility of the 
engine operating in the prescribed manner because the details of the 


burning process are not linked to the engine design and operation. ‘154 


se 


3.2.1.2 Quasi-dimensional Mode] ll! 


Quasi-dimensional models are also structured around a thermodynamic 
analysis of. the engine cylinder during the cycle. Many of the same 
assumptions are utilized to describe the various portions of the process 
as are used in the thermodynamic model. The combustion process, on the 
other hand, is based on more fundamental physical quantities such as 
turbulent intensity, turbulent mixing, jet characteristics in jet mixing 


and the kinetics of the fuel-oxidation process. 
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The quasi-dimensional models can be utilized for the same purposes 
as the zero-dimensional models except that they can now be used where 
changes in the combustion process can be a dominant factor. The major 
drawback of the quasi-dimensional model is its inability to examine, in 


detail, the interaction between fluid flow and engine geometry, (14) 


3.2.1.3 Multi-dimensional Mode 13 


In a multi-dimensional model, the governing partial-differential 
E ions describing conservation of mass, momentum, energy and species, 
and the sub-models describing turbulence, chemical kinetics, and etc. 
are numerically solved subject to boundary conditions and other 
restraints. These models have the potential for examining the 
interaction between fluid flow and engine geometry that is lacking in 
the quasi-dimensional model. The detailed model will predict engine 
performance and emission characteristics from the first principles with 
virtually no empirical relationships. Unfortunately, solving the 
relevant conservation equations in three-dimensional, time dependent 
formulation, coupled with the state equations and sub-models leads to a 


computer program that will tax even the most capable computer system. 
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Chapter 4 
THE TWO CYLINDER CYCLE COMBUSTION MODEL 


Since the two cylinder cycle reciprocating engine is a totally new 
concept, combustion modeling can ke even more difficult than for a 
compression ignition engine. However, the approach taken models the 
expansion cylinder of the new cycle after a diesel engine cylinder. The 
beginning of the expansion stroke will simulate a diesel engine with its 
piston at TDC with a charge of air. For this initial combustion model, 
the air will be assumed to be contained within the cylinder, at 
pressure, with no additional air added after expansion, as in the actual 


new engine cycle. 


4.1 Assumptions 


The assumptions for this single zone combustion model are 
essentially those previously mentioned for the thermodynamic type of 


models. 


a. The First Law of Thermodynamics is used to establish an energy 


balance to determine the temperature at the end of each step. 
b. The working fluid is treated as an ideal gas. 
C. The system contents are homogeneous and of uniform temperature 


and pressure. 
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d. The changes in gas properties due to the rate of change of the 


gas composition are considered to be negligible. 
e. Combustion is treated as a reversible heat release process. 


f. Combustion products are formed in the proportions according to 


the law of perfect combustion. 
g. No dissociation of the products of combustion occurs. 


h. Only four gases are considered to be present and are varied as 


required for perfect combustion. 


i. The incoming air charge is assumed to be pure air plus a 


fraction of the residual gases remaining in the cylinder. 


4.2 Thermodynamics of Internal Combustion Engines 
4.2.1 Ideal Gag {16} 


The assumed thermally ideal gas obeys the state equation 
pV = MRT 
where p = pressure 
V = volume 
M = number of moles 
R = universal gas constant 


and T = temperature. 


The specific gas constant, R, can be written in terms of R and m 
the molecular weight of the gas. 
R= of 


ocr 





If the mass of the gas, m-Mm , then the state equation can be written 


as: 


pV = mkT. 


The specific internal energy for an ideal gas can be represented as 


a function of temperature: 
u= £(T) 
where u = specific internal energy 
and f(T) = function of temperature dependent on 


the gas. 


If the function f(T) is expressed in the form of a limited power series, 


E) 


nas 


alt‏ نا = ا 
where a, to ag are constants which vary‏ 
depending on the gas‏ 


and ü, = internal energy at absolute zero. 


The specific heat at constant volume can be defined as: 
Cy - (dq/ dT) „ =(du/dT) 
Thus, following the same procedures as for the internal 


E 07,18) 
n: 


G = Z pag 
ne 
The specific enthalpy, h, for an ideal gas is given by: 
h= u + Rf. 
It follows that; 17/18) mee: 
= h(T) =u, + 7 AT + ۰ 
ml 


- 29 


energy, 





At absolute zero, T=0: 

h = hy = up. 
Therefore, for a perfect gas, the internal energy varies linearly with 
temperature as: 


no ho DS RT. 


The specific heat at constant pressure, Cor is defined by: 


C = = ° 
p (dq/dT) (dh/ aT) 5 
For a perfect gas: 
Co =C, +R. 
Now, enthalpy can be expressed by: 
h= Ay + Cor. 


For thermodynamic processes with gases of constant composition and 


specific heats undergoing state changes; 


h. Eus 0. 
Then, 
u= CT 
h = CT: 
p! 


h-u= (C5 =- Cy? Y = RT; 


and Ç, - Cy = R. 


Gas data are often given in terms of enthalpy vice internal energy. 
The conventional form is: 


h(T) /RT 


NZ ho) /RT 


4 
ay + a? + ag” + a,” + asT . 
and the internal energy is expressed as: 
= 4 
u(T)/RT = (a, -1) + aT + ajT* + sr + aT’. 
-30- 





The values for the polynomial coefficients, ag to ag are provided in 


Table 2. 


specific heat are available in the literature. 


zl 
Co» 3,0959 


HO 3.4292 


2.13114E-03 


5.65590E-04 


23 


4 .95240E-0 8 


O, 3.25304 6.52350E-04 -1.49524E-07 
N, 3.34435 2.94260E-04 1.95300E-09 
CoH g -0.71993 4.6426 E-02 -1.68385E-05 
CH, 1.13711 1.45532E-02 -2 „95 876 E-06 
4.2.2 Properties of Gas Mixtures ®? 


Mixtures of gases obey the following. 


a. The gas mixture as a whole obeys the equation of state, 


8.66002E-11 
- 1 
121 9۹-1 
-6 .57470E-12 
-2.67009E-09 


0.0 


pV = MRT, where M is the total number of moles of all species. 


Other formulations for the calculation of enthalpy and 
150267275 


13 


0.0 
0.0 
0.0 
0.0 


0.0 


b. The total pressure of the mixture is equal to the sum of the 


pressures which the individual components/species exert. 


c. The internal energy, enthalpy and entropy of the mixture equals 


the sum 


of the internal energies, enthalpies and entropies which each 


individual component/species would have if it separately occupied the 
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entire volume of the mixture at the same temperature. 


Thus, for mixtures of ideal gases the mole fraction is given by: 
where M, = moles of a specie 
and M = total number of moles. 

Then, 

DX. = 1.05% 
Enthalpy is given by: 

H = 2 بط با‎ = M 2 x;h;. 
Internal energy is given by: 

U- 2Mu, - ME x,u,. 


Specific Heats are given by: 


cp ^ * XiCpi 
Cy P EXC 
4.2.3 The First Law of Thermodynamics*??! 


The emphasis of this model is the closed portion of the cycle. 
Therefore, the First Law of Thermodynamic for closed systems is simply: 
| dQ - dW = ٥ن"‎ 
where dQ = heat energy transfer 
dW = work energy transfer 
aU‘ = change in internal energy. 
The internal energy is defined by: 
U' = U + KE + PE 


where U = the intrinsic internal energy 
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KE 


kinetic energy 


PE = potential energy. 


For a closed system, we can assume that PE = KE = 0. Therefore, 
dQ - dW =dU 
where U = M > تا وک‎ 
M = total number of moles 


X. 
L 


mole fraction of gas i 


U. 
1 


specific internal energy of gas i. 


For non-reacting closed systems, we can write: 
dQ - dw = dU 
where dW = pdV = ( <2 x;p)dV 


and dU = MdCExiuj. 


For a reacting closed system, we can expand this to: 
dQ - pdv = (U, 


where Un 


AU, = heat of reaction 


SU UO U 


- Eo! = GER 


P 


U (T) = energy of products as a 


P 
function of time 
U,(T = energy of reactants as a 
function of temperature. 
4.3 Transfer from the gas to th inder 


To be able to balance the energy in a real system, the heat 
transfer from the combustion gas to the walls of the cylinder must be 
considered. Two basic equations are generally accepted for use in cycle 
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calculations. These are the correlations developed by Annand and 


Woschni. The relationship by Woschni (19) is based upon a forced 
convection model. 
y aca? — KN TT) 
where Ci C, and C4 = constants 
A = area 
D = cylinder bore 
p = pressure 
T_ = mean gas temperature 
= wall temperature 
= piston velocity 
peus motoring pressure 
p' » trapped pressure 
V' s trapped volume 
T' = trapped temperature. 
Although Woschn1's expression is readily accepted, it does not 


separately distinguish between convection and radiation. 


The Annand equation is also largely based on turbulent convection. 
Unlike the Woschni correlation, Annand claims that the Reynolds number 
is the major parameter affecting convection. Convection is the first 
term in his equation. The second term in Annand's equation is a 
radiation term assuming grey body radiation. Thus; 120? 

YA = a(k/D) (Re) (TF, + c(Tg*-Tw*) 
where q = heat transfer rate 


A = area 


EE 





a,b,c = constants 

k = thermal conductivity 

D = bore 

Re = Reynolds Number = pVED/u 
p = density 
Vo = piston velocity 
u = viscosity 


temperature of gas (mean) 


Tg 
Tw 


temperature of wall. 


The range of values for Annand's constants are: 
for a four stroke engine: 
a= 0.26 


b = 0.75 + 0.15 


3.88 21.39 x 109  J/sm*K* 


@ 
for a two stroke engine: 

a = 0.26 

b = 0.64 + 0.10 

c = 3.03 £ 1.06 x 10^? 


Since Annand's equation separates the convective term from the 
radiation term, it is believed that the Annand correlation is better 


Suited to the new cycle calculations. 


4.4 The Combustion Model 
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In the process or heat release from combustion, both physical and 
chemical effects are involved. Liquid fuel injected into an engine 
must be heated, vaporized, and mixed with oxygen in the preparation 
process prior to combustion. Once the fuel is prepared, it may then 
burn at a rate controlled by chemical kinetics. It has been 
demonstrated that the time required for combustion of the- prepared fuel 


is negligible as compared to the preparation time. 


At the beginning of the burning period, chemical kinetics are 
important due to the low temperatures. When fuel is first injected into 
a Cylinder of a diesel engine, the temperature is generally such that 
rapid burning will not occur. Additionally, the heat transferred to the 
incoming fuel causes the temperature to drop in the cylinder. As the 
temperature rises in the cylinder, the combustion rate rises, thus 
increasing the temperature. The heat release rate continues to rise 
until the lack of prepared fuel becomes the controlling factor. When 
the excess prepared fuel is depleted, combustion will proceed at the 
rate of fuel preparation. Figure 12 represents the effects of 
preparation rate and reaction rate in premixed burning as a function of 


crank angle. 


4.4.1 Preparation of Fuel 


After injection, the fuel is physically prepared for combustion. 
As mentioned before, this process involves the atomization, vaporization 
and mixing of the fuel with air. The rate of preparation can be assumed 
to be proportional to the total surface area of the fuel spray droplets. 
If all the droplets are assumed to be of identical size, then it 


E 





follows; 7:14:21? 


Mi 


M 


np 7 DZ 6 


np z D 


where M; = Mass of fuel injected 
M, = Mass of fuel unburned 
n = number of fuel droplets 
p = fuel droplet density 
D. = Initial droplet diameter 
D = Droplet diameter. 
The total area 
Area = nxD2 = n z (6M /np=< ر‎ 
Area = (6M,/npD.2)1/3 (6M /pz)2/3 
Area = E Zuep 


Assuming that the density, p, and initial diameter, Dor 


are constant, 
then the 


Area a M.1/3y 2/3, 
1 u 


- 


Allowing for the effect of oxygen availability on the mixing of the 
fuel, the preparation rate, PR, can be written as: 
_ l-X, x5 m 
PR z KM; M, Poo 
where x = empirical constant 


m = empirical constant 
Poo = partial pressure of oxygen 
K = constant. 
The constant K is a function of the characteristics of fuel injection, 


air movement and combustion chamber shape. Typical values for four 


n 





stroke engine are; 143 


K = 0.008 - 0.020 
x= 2/3 


m= 0.4. 


4.4.2 Reaction of Fuel 


Since diesel fuel is not a pure substance, it is impossible to 
ascertain the exact chemical equations involved since the actual 
compounds in the fuel are unknown. The temperatures that are available 
from experiments are only average cylinder temperatures. With these 
approximations/estimations, the equations for reaction rate are highly 
empirical. The degree of approximation involved may be justified due to 
the short time period during which chemical kinetics is of importance. 
Also, the total fuel that is burned is equal to the amount of fuel that 
is prepared. The reaction rate equation that was proposed by Whitehouse 


(7,14,21) is based on the Arrhenius equation. 


and Way 
R = (K'Pog)/ (MVE ) | (PR-R ax exp (-act/T) 
where R = reaction rate per degree crank angle 
K' = empirical constant 


act 


empirical constant 
Poo = partial pressure of oxygen 
PR = preparation rate 


N 


engine speed in rpm 


T 


cylinder temperature, 
The effect of the ignition delay period is incorporated in the Arrhenius 


type expression exp(-act/T). Typical values of K' and act are : 
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act = 1.4 x 104 
K` = 1.2 x 1010 for two stroke engines 
K` = 65 x 1010 gor four stroke engines 


4.5 Verification of the Model 


The model was converted top computer code using TRS-80 Model III 


Disk Basic. The program listing is presented in Appendix B. 


In an effort to set the empirical coefficients, the average value 
was used for all coefficients that had a range of values for four stroke 
engines. The program was run and compared to the data obtained by 
Remley +22? in actual engine testing in the Sloan Automotive Laboratory. 
Figure 13 represents the pressure versus volume curve for the model and 
for the engine run by Remley. Appendix A provides specifications of the 


test engine. 





Chapter 5 
SELECTION CF FUEL INJECTION POINT 


In order to obtain the maximum work and highest efficiency from the 
new two cylinder cycle, the time of fuel injection should be optimized. 


To obtain this optimum, a number of cycles were run on the computer. 
5.1 Selection of the Mod fficients 


The model was run assuming the expansion cylinder at TDC with an 
air charge at a temperature and pressure of 1090°K and 10 atmospheres 
while the engine speed of 850 rpm and air/fuel ratio were held constant. 
The selected fuel was CoH} g (iso-octare) with a lower heating value of 


4.2 x 10’ joules/kilogram and a residual air fraction of 0.05. 


The model was run several times to obtain a value of K in the 
equation: 
PR = rm, (IX) m Xp. 
The values of x and m were held constant at 2/3 and 0.04, respectively, 
as the values used for four-stroke diesel engines. When searching for a 
value of K, a diffusion combustion period of 70 - 120 degrees of crank 
angle was sought. This was found through several iterations to occur at 


a value of K = 0.012. 


ES 





The values of constants for the reaction rate equation: 
R= [K'Pog/NT? "Š lexp(-act/T) | (PR-R) dx, 


were selected as the values for four-stroke diesel engines, 


With this input data and selection of constants, the model yields a 
heat release rate curve which closely resembles that described by 
Ricardo, Lyn and Whitehouse et al, see figure 14. The premixed burning 
phase yields approximately 45 percent of the heat release, the diffusion 
controlled burning phase yields approximately 45 percent of the heat 


release with the tail of combustion providing the remaining 10 percent. 


5.2 Optimizing Fuel Injection 


Intuitively, the maximum work and highest efficiency would be 
expected with fuel injection and combustion occuring at TOC, or 
immediately thereafter. This, however, does not appear to be the case 
when the data is evaluated. See figures 15 through 21. While the fuel 
injection is varied from 180 (TDC) to 205 degrees crank angle, with an 
injection period of 20 degrees, the thermal efficiency rises. For fuel 
injection occurring from 180 to 195 degrees, the temperature at 360 
degrees (BDC) is not sufficient to heat the regenerator matrix to a 
temperature which will pre-heat the incoming air charge to 1090 degrees 
Kelvin as specified by the input data. For fuel injection occuring at 


205 degrees, and later, incomplete combustion will result. 


From this approach, the optimum point of fuel injection occurs at 


200 degrees crank angle for a fuel injection period of 20 degrees. 
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When the fuel injection period is reduced to 10 degrees, a similar 
pattern is observed. A fuel injection point on, or before, 200 degrees 
results in the cylinder gas temperature dropping too low to support 
sufficient air charge pre-heat. Fuel injection on, or after, 210 
degrees results in incomplete combustion. See figures 22 through 24. 
In this case, the optimum point of fuel injection occurs at 205 degrees. 
The thermal efficiency for this case is higher that the case of a 20 
degree injection period. Also, the specific fuel consumption is lower 


in the case of 10 degree injection as compared to 20 degree injection. 


Through a similar analysis, the case of an air/fuel ratio of 25 
yields an optimum fuel injection point of 195 degrees crank angle for a 
period of 20 degrees. For this air/fuel ratio, the value of K in the 
preparation rate equation was selected as 0.018 to achieve a similar 


heat release rate pattern. 
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Chapter 6 
COMMENTS AND RECOMMENDATIONS 


As can be readily seen, the output from this type of thermodynamic 
model is dependent on the value of the empirical coefficients. The 
characteristics of the heat release rate curve will shift as a function 
of air/fuel ratio, temperature, pressure and engine speed. Therefore, 
only a comparison of results from a defined heat release should be used 


for qualitative comparison analysis. 


Since this computer model was written for a personal computer, the 
time required for one run is excessively long for a detailed comparison 
analysis. The run time for one run with five degree increments is 
approximately 1 hour and 20 minutes. A motoring analysis (no 
combustion) requires approximately 15 minutes. The amount of time 
required in the combustion iteration process is the difference between 
the two. These times were obtained when running the program with no 
remark statements and elimination of all unnecessary spaces in the 
program. Undoubtedly, the efficiency of the program can be somewhat 
increased by utilizing some clever programming techniques. However, the 
use Of a small computer strictly dedicated to a comparison analysis with 
Crank angle increments of one or two degrees can occupy the machine for 


an inordinate period of time. 


AO 





6.1 Recommendations 


This single zone model allows for cycle studies, However, a 
problem that must be explored is the formation of soot and gaseous 
pollutants. This can be accomplished by expanding the model to a two or 


1,123,243 During this model expansion, the effect of 


four zone mode 
chemical kinetics should be further examined to display a more realistic 
combustion process. The values of the coefficients for the polynomial 
expression of enthalpy for the other products of combustion are readily 


available. 117,183 


The effects of heat transfer from the system may be more 
appropriately modeled by the use of the widely accepted Woschni 
correlations, 119? The use of Annand's correlation, however, does allcw 


for the separation of convection and radiation. 


The effects of mixing of the air charge with the fuel must be 
further explored to determine the effects on combustion intensity and 


efficiency. £25} 


The use of a larger computer system would be most beneficial in a 
comparison analysis. Single runs can be easily done on a personal 
computer system, however, many runs using small crank angle increments 
are best, although more costly, performed on a main frame system capable 


of performing numerous simultaneous calculations. 


Lastly, to obtain realistic coefficients for the empirical 
constants in the preparation and reaction rate equations, experiments 
using a rapid compression machine are considered appropriate. This 
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would provide for realistic data with minimum cost. 
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Figure 1: Pressure vs. Volume for New Cvcle vs. Diesel Engine 
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Figure 14: Heat Release Rate Obtained from Computer 
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5٦6 ENGINE COMBUSTION CYCLE 


RUN: CARMICHAEE ENGINE 
INPUT DATA: 
CYEENDeK BORE = S-259 METERS 
STEORE = 25/72: nen 
CONNECTING ROD LENGTH = 745 M 
ENEINE SPEED = 950 RPM 
FNGIHSWRGONPRESSIGN RATIO = 8 S 
۰۲۰۰۸ ۷۳۷۰۸۵٣۲ 0 30 
TRAPPED PRESSURE E 
TRAPPED TEMPERATURE = 1090 DEG KELVIN ` 
RESIDUAL AIR FRACTION = 105 
BUFRZSEEFETEBZEOR THTS ANALYSIS On 122° 7 >20 OCTANE) 
IITA LOVER HEATINC VALUE = 237 2E>07 JOULES/KG 
SEHR AIR / FUEL RATIO = 15.1151 
meee 7 AIR EQUIVALENCE (PHI) = „503855 
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Appendix A 


Specifications of Test Remley Engine “22? 


Type of Engine 

Bore 

Stroke 

Cylinder Displacement 
Connecting Rod Length 
Compression Ratio 

Number of Compression Rings 
Number of Oil Rings 

Number of Inlet Valves 

Number of Exhaust Valves 

Valve Diameter 

Valve Lift 

Inlet Valve Timing open/close 
Exhaust Valve Timing open/close 
Diameter of Intake Manifold Pipe 


Diameter of Exhaust Manifold Pipe 


Four Stroke 
4.0 inches 
2.5 inches 
31.41 cubic inches 


6.25 inches 


14.3 : 1 

2 

1 

2 

2 

1.286 inches 
0.280 inches 
159BTDC/ 509ABDC 
50 BBDC/15%ATDC 


2.00 inches 


1.60 inches 





Data Collected fram Test Run'22! 


Engine Speed in RPM 1450 

Inlet Pressure 13.5 inches Hg gage 
Exhaust Pressure 13.4 inches Hg gage 
Inlet Air Temperature 1869 F 

Air to Fuei Ratio 25.38 

IMEP 88.1 psi 

Start of Injection 12.59 BTDC 

Ignition Delay 5.59 

Period of Fuel Injection Ines 
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Appendix B 
Computer Model 


The computer program is written in TRS-@ Model III Disk Basic and 
consists of a main program and nine subroutines. The program listing 
has numerous remarks statements inserted to make the algorithm and 
computer code easier to understand. Since the program takes a 
considerable length of time to run, it is recommended that the remark 
statements be deleted before running. Samples of output are presented 


in figures 15 to 24, 


== 
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